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In the canonical picture ol the evolution ol the quark-gluon plasma during a high-energy heavy- 
ion collision, quarks are produced in two waves. The first is during the first fm/c of the collision, 
when gluons thermalize into the quark-gluon plasma (QGP). After a roughly isentropic expansion 
that roughly conserves the number of quarks, a second wave ensues at hadronization, 5-10 fm/c 
into the collision. Since entropy conservation requires the number of quasi-particles to stay roughly 
equal, and since each hadron contains at least two quarks, the majority of quark production occurs 
at this later time. For each quark produced in a heavy-ion collision, an anti-quark of the same 
flavor is created at the same point in space-time. Charge balance functions identify, on a statistical 
basis, the location of balancing charges for a given hadron, and given the picture above one expects 
the distribution in relative rapidity of balancing charges to be characterized by two scales. After 
first demonstrating how charge balance functions can be defined using any pair of hadronic states, 
it will be shown how one can identify and study both processes of quark production. Balance 
function observables will also be shown to be sensitive to the charge-charge correlation function in 
the QGP. By considering balance functions of several hadronic species, and by performing illustrative 
calculations, this class of measurement appears to hold the prospect of providing the field's most 
quantitative insight into the chemical evolution of the QGP. 
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INTRODUCTION AND THEORY 



In a central heavy ion collision at RHIC (Relativistic Heavy Ion Collider) or at the LHC (Large Hadron Collider), 
several thousand hadrons are created from the initial collision of a few hundred incoming nucleons. In the central 
unit of rapidity, aside from a few dozen extra baryons, the roughly one thousand hadrons are created and evolve from 
quark-antiquark creation processes. For every up, down or strange quark observed in the final state, one usually finds 
one extra anti-up, anti-down or anti-strange antiquark within roughly a unit of rapidity. For quark-antiquark pairs 
created early in the collision, the balancing pairs might be pulled apart by the initial tunneling process by a fraction of 
a fm in distance, and are then pulled further apart by collective longitudinal flow and diffusion. If the quarks are pulled 
a half fm apart at a time 1.0 fm/c, collective flow would pull them apart by 7.5 fm by the time breakup occurs (~ 15 
fm/c), and diffusion would spread them apart even further. In the canonical view of the quark-gluon plasma (QGP), 
a first wave of quark production occurs when the quark-gluon plasma is created during the first fm/c. The number is 
then roughly conserved in an semi-isentropic expansion until hadronization, when a second wave of production ensues. 
Due to entropy conservation, a thousand partons would be expected to convert to roughly a thousand hadrons, and 
since each hadron has multiple quarks, and since the gluonic entropy has no quarks, the number of quarks should 
more than double during hadronization. If hadronization were to occur later in the process, perhaps at 5-10 fm/c, 
the balancing quark anti-quark pairs created at hadronization would be unlikely to separate by more than a few fm 
before breakup. 

Charge balance functions were proposed as a means for identifying and quantifying the separation of balancing 
charges [1] . They represent the conditional probability of observing a balancing charge in bin p\ given the observation 
of a charge in bin p 2 , and are defined by: 

n / I \ _ (l n +iPi) - n-(pi)][n-(p 2 ) - n + (p 2 )\) m 

B+-[pi \p 2 ) = i , (1) 

{n + (p2) + n_(p 2 )) 

where (n + /-(Pi)) is the probability density for observing a positive/negative particle in bin pi, and (n+(pi)n_ (p 2 )} 
is the probability density for observing a positive particle in bin pi and a negative particle in p 2 . If the number of 
positives and negatives are equal, and if the detector was perfectly efficient for all p\ , integrating the balance function 
over all p\ would give unity. The label pi can refer to any measure of momentum, including rapidity or pseudo-rapidity. 
The observable can be modified to more appropriately treat the case where the net charge significantly differs from zero 
[2]. In short, the balance function is simply the application of a like-sign subtraction with the purpose of statistically 
isolating the opposite balancing charge. 

More generally, balance functions can be analyzed for any two set of hadrons and antihadrons, 

B a p{pi \p 2 ) = -, — r — , (2) 

For instance, one could consider balance functions where a were protons and (3 were negative kaons. The antiparticles 
are noted by a and /3. For this study we will confine ourselves to the situation where the net charges are zero, which 
is certainly a good assumption at LHC energies. For the case where the net charges are not equal, one might wish 
to follow the example in [2j and define the denominator using the lesser of the two charges, n„ or n a , followed by a 
mixed-event subtraction. 

Balance functions can be analyzed in six dimensions as a function of p\ and p 2l though statistics make that 
prospect unlikely. Instead, the condition p 2 is usually the observation of a particle anywhere in the detector, while pi 
refers to the observation of the second particle with relative rapidity Ay, or relative azimuthal angle A0, or relative 
invariant momentum <5i nv . Balance functions were used at the CERN ISR to study hadronization dynamics in pp and 
e + e~ collisions in the 1980s [3HZ]i while their use in heavy ion collisions was motivated by the desire to distinguish 
between early vs. late production of charges pQ. For more central collisions, balance functions have been observed to 
significantly narrow when binned in relative rapidity [8 10 . This behavior is quantitatively consistent with the idea 
that a good fraction of the charge is created late in the collision, as expected from delayed hadronization with the 
existence of a long-lived quark-gluon plasma |llj . Narrowing is also predicted and observed as a function of relative 
azimuthal angle [TT1 [TZ] , though this paper will focus on the behavior in relative rapidity. 

Our first goal is to understand how to calculate balance functions between any two hadronic species a and /3. 
We especially wish to know what happens if charge production comes in two waves, the first wave being the initial 
thermalization of the QGP, which is where the quark number rises quickly from zero in the first ~ 1 fm/c, and the 
second wave at hadronization, which may be in the 5-10 fm/c window. This goal is complicated by the fact that 
hadrons carry three charges (one for each light flavor of quark or equivalently strangeness, electric charge and baryon 
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number) , which makes the problem rather entangled. Charge balance functions binned in relative spatial rapidity are 
characterized by a scale ct(qgp) before hadronization, which might well be the greater part of a unit of rapidity. During 
hadronization a second group of balancing charges is created with relative coordinates characterized by cr^ad) ~ 0. In 
the next section, it will be shown how one can use local charge conservation overlaid onto an assumption that extra 
charge within a volume is distributed thermally to derive expressions for both components of the balance functions. 
From this perspective, all balance functions in relative spatial rapidity will be determined in terms of cr^d) an d 
ct(qgp), the number of quark species per unit rapidity before hadronization, dN u ^ s /dy, and the number of hadrons 
per unit rapidity, dN a /dy in the final state. The hadronic yields are experimentally measured and the rapidity density 
of quarks can be estimated from the total entropy, thus leaving the widths as the least understood quantities. The 
balance function in terms of the relative spatial coordinate along the beam axis translates into a balance function 
in relative rapidity after convoluting with a thermal kinetic distribution and including decays. In the final section, 
we present illustrative predictions for balance functions in relative rapidity for several species by using a thermal 
blast-wave model to map between spatial rapidity and momentum space rapidity. 



THEORY 

Balance functions are related to charge correlations. For the purposes of this derivation the charge densities are 
considered as a function of the coordinate 77, which describes the longitudinal position in Bjorken Coordinates. 

z = rsinh(?7), t — TCOs(ri) 7 (3) 
t = \J t 2 — z 2 , r] = tanh - 1 (z/r). 

In the absence of longitudinal acceleration a particle moving with the fluid has fixed 77, and aside from diffusion the 
separation of balancing charges would be fixed in A77. 

Before progressing, we define the charge correlation function, 

9ab{m,m) = (p a (m)pb(m)y (4) 

where the sum over i ^ j covers all particles i and j, and the prime notes that the correlation of a charge with itself 
is subtracted. The indices a and b refer to the specific charge, e.g. net strangeness. For this paper Roman indices will 
refer to charges, e.g., the net number of up, down or strange quarks, while Greek indices denote specific species, e.g., 
7r + ,p, K~ . Since a chargeless plasma is being considered, (p a ) = 0, and one need not subtract the terms (p a }(pb)- For 
a hadronic system conserving baryon number, electric charge and strangeness, the index can equivalently sum over 
the net number of up, down and strange quarks. The charge-charge correlation can be expressed as: 

9++(vi,V2) = {[n+ivi) - n-(»7i)][n + (7ft) - n_(r? 2 )]}' (5) 
= -B+-(vi\m) (n+(m) + n-(V2)) ■ 

Here, the positive and negative subscripts refer to the sets of all positive, or all negative particles. The relation 
between the balance function, B a p and g a b becomes complicated if particles have more varied charges, which is the 
case for a hadronic system, e.g., the E~ carries baryon number, electric charge and strangeness. These cases will be 
discussed in the next few paragraphs. 

The reason we switch from balance functions to correlations is that the correlation does not change suddenly at 
hadronization, except for where rji — r\i- This follows from local charge conservation. It can be understood by 
considering the addition of a pair with r\\ ~ 772 during hadronization. The contribution of this single pair to the sum 
in Eq. ^ where either i or j points to any other particle besides those from the created pair vanishes, because one 
is considering the creation of a pair with equal but opposite charges at the same point. The only contribution comes 
from the element of the sum where both particles come from the pair, which then shows up at 771 =772. Assuming 
hadronization is sudden, and assuming one understands the charge correlation before hadronization, one would also 
know 3(771,772) immediately after hadronization, except for the region 771 w 772. To determine the correlation in the 
region of small relative coordinate, one can use the sum-rule for charge correlations, which follows from integrating 
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the definition of g in Eq. Q , 

- J dr)igab(r)i>Vi) = ~ Y2 Qj,aS(V2 - Vj)qi,b (6) 
j 

= Xab{m) = ^2(n a (v2))q a Mq a ,b- 

a 

The first step used charge conservation, J^. % = 0. The average number of particles of a given species a within dr\ is 
(ji a {r\j)dr\. Assuming instantaneous hadronization, in order to satisfy the sum rule of Eq. ([6]), the charge correlation 
immediately after hadronization must be: 

9ab{m,m) = 9c2 GF \m,m) + 9^ ivum), ( 8 ) 
g { ab d) {m,m) = - [x { aT d \m) - xi? GP) (m) <5(m - m), 



(had), 
tab \ 

,(QGP), 



\„;, (V)= q^aq a ,b( n a (v)) 

\„/; (v) = qa,aq a ,b(n a (v)) , 

QSQGP 



where g(Q GP ) describes the correlations both immediately before and immediately after hadronization, but neglects 
the hadronization component created at 771 — r\2- The sums over a cover the species for each state, i.e., over partonic 
species for the QGP state and over hadronic species for the hadronic state. The value \ab, when multiplied by the 
delta function, represents the charge-charge correlation that would ensue from independent particles, i.e., when the 
only correlations come from a particle with itself. Here the values (n a ) are the densities per unit 77 of the species a, so 
if one measures the final-state yields xift"^ can be considered as known. The values of Xab can also be extracted from 
a one-body treatment such as hydrodynamics. The matrix Xdjj is diagonal in a QGP if the charges refer to the 

net number of up, down and strange quarks. In contrast, hadrons have multiple charges and xl& ad ' ) nas off-diagonal 
elements. Since hadronization is sudden, but not instantaneous, one would expect to replace the delta function with 
some function of finite but narrow width, normalized to unity. 

Our next goal is to determine the balance function for any hadronic species just after hadronization, given g a \, in the 
QGP phase. Eq. Q describes how to extract g a b just after hadronization. However, once there are multiple charges 
spread across a variety of species, it is not easy to understand how the correlation functions, g a b{Vi-i r ]2), determine 
the balance functions, B a p(rii\rj2). Here, a and b refer to any conserved charges, while a and /3 refer to the charge 
carried by a specific species, where the particle and anti-particles of each species are denoted by a and a, 

B a 0{m\m) = / — , , , — 7—T\ ( 9 ) 



Here, n a is the density (number per unit rj) of particles of species a. Thus, g a /3 is the correlation of the effective 
charge defined by the number of a specific species minus the number of its antiparticle. With this definition, one can 
see that 

9aP = ~9a,/3, 9&Si = -9afi, 9a^=9afi, (10) 

B a p = ~B a p, Ba.fi = -B aj 3, B a = B a p. 

As an example, one can consider the proton- K ~ balance function. In this example, the index a would refer to protons 
and /3 would refer to negative kaons. The corresponding charge correlation function would be 

9 P K-{m,m) = (KO71) - n p (r]i)][n K -(r]2) - n K +(r]2)]). (11) 

The suffixes a and /3 can also refer to a subset of species, with a and j3 referring to the equivalent subset of 
antiparticles. For instance, a could refer to the set of all positive particles, while f3 could refer to the set of all 
antiparticles. Switching the indices leads to the relations: 

9af3 = 90a, BapUp = Bp a n a . (12) 
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Determining the balance functions for arbitrary species requires making the jump from g a b to g a p. There are three 
conserved charges, which we will consider to be the net numbers of up, down and strange quarks. Although one could 
have equivalently used baryon number, electric charge and strangeness, the quark numbers are more convenient since 
one does not expect any off-diagonal elements to g a b in this basis for the QGP. For the species-labeled correlations, 
ga/3, there are many more possibilities in the hadronic state. Even for the final state, one might wish to consider 
charged pions, charged kaons, protons or lambdas. Neutral kaons must also be taken into account for absorbing 
strangeness, but because they oscillate into K s and Ki, cannot be easily used for balance functions. Since g a p has 
more elements than g a b, additional assumptions are required if g a p is to be determined from g ab . 

Observing a hadronic species a at position 771 infers one has observed the three charges q a>a , which is the number 
of up, down and strange quarks in the resonance a. The correlation gabivii^) should then provide the probability of 
finding the balancing charges at position rj 2 . In order to deterine g a p one then needs a model to determine how an 
extra charge qb at position 772 influences the probability of finding a hadronic species j3 at the same position. 

By assuming that the local distribution of hadrons is determined by a thermal distribution constrained by the 
local charge density, one can determine g a /3 from g a b- To show this, we express the two particle correlation as being 
determined by a grand canonical ensemble with Lagrange multipliers applied to constrain reproduction of the average 
two-particle correlation function, i.e., 



(AB) = ilJ ABe-f d " H «/ T ^exp f d m d m Y,PMHab{r}i,m)pb{m) 

I L ab 

dr]idr]2 ^2 PaM^abim, m)Pb{m) 



(13) 



Z = TJ e -^" H °/ T Wexp 



Here, Hq is the Hamiltonian or relevant free energy density, T is the temperature, and Ha^iviiVz) plays the role of a 
Lagrange multiplier chosen to enforce that g &(?7i7 ^2) is reproduced. The strategy will be first to find p a b in terms of 
g a b, then to use p a b to determine g a p. The correlation function <7 a b(??i, 772) is found by replacing the operators A and 
B above with 



a p 



(14) 



where a and (3 are summed over all hadronic species. By assuming that the weighting is proportional to an exponential 
of the constraint (fixing g a b), this is essentially a thermal ansatz. 

Since the correlation would be zero if not for p, we can expand the expression for small p and find: 



9ab(m>V2) = ^2{n a (rii))q a , a qp,b{np{r)2)) exp < ^2q a ,cPcd(m,V2)qp,, 



(15) 



a/3 



cd 



(n a (i]i))qa,aqa,cPcd(m, r l2)q0,dqi3,b(ni3{ri2)), 



ct(3cd 



^2 Xac{m)pcd(m,V2)xdb{m), 



cd 



where \ was defined in Eq. ([6]). The assumption of small p is warranted given that charge-conservation correlations 
are small (at least for central collisions). Inverting the equation, one can then find p a b in terms of g a b, 



pab(m,m) = ^2xL 1) (vi)9cd{m,'n2)Xdb (%)• 



(16) 



cd 



One can now find g a p by inserting 

A = n a (r)i) - n s (r)i), B = npfa) ~ ^(772), 



(17) 



into Eq. (13 1. Here n a is the density of the anti-particles to a. Again, assuming equal numbers of particles and 
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antiparticles, (n a ) = (n a ), and assuming that [i ao is small, 

g a p(vi,V2) = ([n a {m) ~ na(Vi)} [M^a) - "§(772)] ) (18) 
= (n a (m)) (npim)) cxp < ^2q a MfJ-ab(m,m)lp,b > + (n a (?7i))(n^(77 2 )) exp < ^2q a , a fJ, a b(r}iim)qp,b > 

{ ab ) { ab J 

-(n a {r]i)}( n p(V2)) exp I 

I ab J I ab J 

~ 4(n Q (ryi))q Q ^^ Qb (77i, l) 2 )9fti(^())2)). 
From Eq. ([9]), one then finds an expression for the balance function, 

£0^(7711772) = 2 ^(n a (77i))g a!a /x a 6(77i, 772)^,6 (19) 
= 2^(n Q (77 1 ))g ct!a xi7 1) (?7i)ffcd('7i,?72)xib 1) ( 7 72)9/3,b- 



abed 



One test of this result is to see whether integrating the balance function over all 771, summing over all a, and weighting 
with q a , a , one should get the net amount of charge a found in other particles due to the condition of having observed 



a particle of species (3 at position 772. Performing these operations from the expression for B in Eq. (19), 



^2 / dr ii q a ,aB al5 {-qi\ri2) = 2 j dr^ ^ q a , a {n a (r]i))q atb xi c 1 \m)9cd(vi,V2)x d b (V2)q^,b (20) 

a J J abed 

= 2 dru Xab{vi)xic~ 1) (vi)gcd(m,m)x ( db 1) ( f i2)qf3,b 



dm gac(m,V2)xib 1) (v2)qfs : b 

-2 Xac{m)x ( cd 1] (V2)qp,d 

cd 

-2qp, a - 



The second-to-last step used the sum rule for integrating g in Eq. pi). The factor of two comes from the fact that 
the sum over all species, a, double-counted the contributions. For instance, the term for which a — tt + also includes 
the contribution from ir~ , and the term for a = tt~ also includes the contribution from the ir + . 

CALCULATING WEIGHTS FOR BOTH COMPONENTS FOR ALL HADRONIC SPECIES 

From Eq. ([8]), one expects two components to the charge correlation <7ab(7?i, 772). Assuming a boost- invariant system, 
one can assume a dependence on A?7 = 77! — 772 , rather than on 77! and 7/2 individually. This expectation inspires one 
to write the balance function for all species B a p(Arj) in terms of two components, 

B a p(A V ) = w^ GP) b^ GP \A v ) + w^¥^(A V ), (21) 

where 6^ GP - ) and &( had ' are both normalized so that J dArjb(Arj) = 1. 

The weig hts, W (Q GP ) and «;( had ), can be determined from the charge correlations, which in turn depend on the 
matrices \ab- From Eq. (fsl) , 



fj„ 



ab (A v ) = x ^ GP) b^(A v ) + M h fc ad) - xi? GP) l b^(^). (22) 



Here, the delta function in Eq. ([8]) was replaced by a Gaussian of finite width, where the width is determined by 
the charge diffusion between hadronization and breakup. The correlation before hadronization, g^ , should be 
diagonal if quarks are good quasi-particles, 

xi?° P) = (n a +n- a )6 ab , (23) 
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where n a is the density of up, down or strange quarks, and n a is the density of the antiquarks. In this formulation 
there is an explicit assumption that the diffusive widths of the charge correlation before hadronization are independent 
of flavor. Whereas the form of ^(Q GP ) is model dependent, = ( n a)Qa,aQ0,b is determined from final-state yields. 



in Eq. (21) 



After inserting the above expression for g ab into Eq. (19), one obtains B a p, from which one can read off the weights 



- oSTlr, \n — 1 (had) (QGP) — 1 (had) 



(QGP) _ \ -l(had) (QGP) -l(had) 

W a/3 

abed 

(had) _ 9 V^/ v -l(had) ,„(QGP) 

ab 



The characteristic width of &(Q GP ) is determined by the charge correlation before hadronization, g^ GP \Ar]), and 
one might expect it to be of the order > 0.5. In contrast, &( had ) is characterized by a narrow width describing the 
diffusion of charge after hadronization and might have a width ~ 0.1 — 0.2. Although the derivations assumed that the 
species were locally populated according to local thermal equilibrium, the weights are completely determined given 
the populations for quarks just before hadronization, and the rapidity density for hadronic species (n a ). 

Whereas the hadronic populations can be taken from experiment (or from a thermal model tuned to experiment), 
the number density of quarks just before hadronization is dependent on model assumptions. Even if one uses en- 
tropy arguments to infer the number of quarks, neglecting the entropy created during hadronization, the number of 
quarks can depend on how much entropy was carried by gluons. For that reason the ratio of the rapidity density 
of quarks before hadronization to the rapidity density of final state hadrons was varied. Three ratios were explored: 
^quarks /^had =0.7, 0.85 and 1.0. The hadron density included neutral hadrons, and the decay products of strange 
baryons and the K s . 

Despite the wide coverage and detailed analysis of RHIC data, the uncertainty in the yields of particular species 
at RHIC can be rather large. Whereas the yields of pions are known to better than the 10% level, yields of protons 
and anti-protons are uncertain at the 25% level. Given these uncertainties, we use yields from a thermal calculation 
based on a temperature of 165 MeV. The calculation involved generating particles thermally from a hydrodynamic 
evolution. Particles of all hadronic species were then evolved through a hadronic cascade, whose main purpose was 
to model the hadronic decays. Since only the yields were sought, the dynamical evolution of the cascade was rather 
inconsequential. Weak decays were not performed. The remaining species and their yields for central collisions are 
given m Table [TJ The yields given in the table were then modified by an additional factor fg , which reduced the yields 
of all baryons by the same factor. Given that the number of anti-baryons is less than the number of baryons at RHIC 
by a factor of 0.7, and given that only anti-baryons are accompanied by an additional charge, one might expect to a 
factor of fs ~ 0.85. Comparing the numbers below to proton yields from PHENIX [TJ], one would expect fg as 0.5, 
whereas a value closer to 0.7 might be expected from STAR's yields [UJ [TS]. The cascade code did not include 
baryon-baryon annihilations, and it is not clear whether a thermal calculation followed by hadronic processes would 
be more consistent with the STAR or the PHENIX values. Given these uncertainties, calculations were performed for 
three values, fs = 0.5, 0.6 and 0.7. 

Only a partial list of resonances were applied in the calculation of weights, those that survive to the final state aside 
from weak decays. For baryons this includes protons, neutrons, £ + /~, EE - / , Q~ and the corresponding antibaryons. 
The included mesons are 7r + /°/ - and K + '~'° . The thermal model provides yields, (n) a , at midrapidity, and are listed 
in Table [I] The "thermal" model was a hydrodynamic model followed by a cascade simulation, where hadronization 
was performed thermally with a temperature of 165 MeV. 

The resulting weights for the default calculation (f B = 0.6, ^ qU arks/«had = 0.85) are shown in Table |h} The 
weights for the tt + tt~ balance functions were not surprising. In the default calculation the hadronization process is 
responsible for nearly two thirds of the final quarks (2 for a meson plus three for a baryon). Pions represent ~ 80% 
of the final-state particles, and given these facts it was not surprising that the hadronization component of the 7T7t 
balance functions integrated to 0.636, while the QGP component integrated to 0.239. The sum did not integrate to 
unity because observing a charged pion does not ensure that the remainder of the system has one extra pion of the 
opposite charge due to the possibility of the charge being balanced by other species. If one could measure the balance 
function in coordinate space, i.e. as a function of rj, there would be a large narrow peak from the hadronization 
component and a smaller broader structure from the QGP component. 

The default results for the K + K~ balance functions were also in line with expectations. Since rather few addi- 
tional strange quarks are produced during hadronization, the hadronization component turned out to be quite small. 
Observing the lack of a narrow peak in K + K~ balance functions would confirm the notion that the QGP was indeed 
rich in strangeness. 



hadron species 


yields, n a 


p,n,p,n 


22.5 


A, A 


8.5 




3.4 


n- c° S-o 0+ 


1.95 




0.35 


+ - 
7T , 7T , 7T 


268 




54 



TABLE I. Hadronic yields used to calculate weights. Yields were calculated from a hydrodynamic/cascade model with no net 
baryon number, where the initial hadronic populations were set according to a thermal distribution with a temperature of 165 
MeV. An additional factor, /s, was applied to the baryon yields listed here to account for the experimental uncertainties and 
for greater consistency with experimental observations. 





P 


A 


E+ 


£" 










K+ 


V 


0.441,-0.066 


0.485,-0.162 


0.491,-0.146 


0.479,-0.178 


0.535,-0.242 


0.529,-0.258 


0.578,-0.338 


0.006, 0.016 


-0.044, 0.096 


A 


0.183,-0.061 


0.242,-0.094 


0.242,-0.094 


0.242,-0.094 


0.302,-0.128 


0.302,-0.128 


0.361,-0.161 


0.000,-0.000 


-0.059, 0.033 


fr 


0.074,-0.022 


0.097,-0.038 


0.099,-0.033 


0.095,-0.043 


0.122,-0.049 


0.120,-0.054 


0.144,-0.064 


0.002, 0.005 


-0.023, 0.016 


E+ 


0.072,-0.027 


0.097,-0.038 


0.095,-0.043 


0.099,-0.033 


0.120,-0.054 


0.122,-0.049 


0.144,-0.064 


-0.002,-0.005 


-0.025, 0.011 




0.046,-0.021 


0.069,-0.029 


0.070,-0.028 


0.069,-0.031 


0.093,-0.036 


0.092,-0.038 


0.115,-0.045 


0.001, 0.001 


-0.023, 0.008 




0.046,-0.022 


0.069,-0.029 


0.069,-0.031 


0.070,-0.028 


0.092,-0.038 


0.093,-0.036 


0.115,-0.045 


-0.001,-0.001 


-0.023, 0.007 


n+ 


0.009,-0.005 


0.015,-0.007 


0.015,-0.007 


0.015,-0.007 


0.021,-0.008 


0.021,-0.008 


0.027,-0.009 


-0.000,-0.000 


-0.006, 0.001 




0.119, 0.318 


0.000,-0.000 


0.239, 0.636 


-0.239,-0.636 


0.119, 0.318 


-0.119,-0.318 


-0.000,-0.000 


0.239, 0.636 


0.119, 0.318 


K~ 


-0.175, 0.384 


-0.627, 0.352 


-0.603, 0.417 


-0.651, 0.288 


-1.055, 0.385 


-1.079, 0.321 


-1.507, 0.354 


0.024, 0.064 


0.452, 0.031 



TABLE II. Default results for the weights w^i , resulting from the thermal model described in the text. The weights 

describe the contribution to the balance functions B a p(Ari) from the correlations driven by the charge correlations just before 
hadronization and the additional correlation that appears during hadronization and is local, Arj ~ 0. 



The pp balance function came out contrary to expectations expressed in previous papers [T] . Even though protons 
are composed entirely of up and down quarks, and even though a large fraction of up and down quarks are produced at 
hadronization, the hadronization component is small, or perhaps negative. This comes from the fact that the strength 
of the hadronization component, as determined by the sum rule in Eq. depends on the density of observed 

baryons. Since the observed number of baryons is rather small, the sum rule can be saturated by the number of 
baryons in the QGP component. If one were to consider baryons alone, the sign of the hadronization component in 
the baryon-baryon correlation depends on the sign of Xu, ~ xi? GP \ where bb would refer to the baryon charge. 
Since the baryon number of a single quark is 1/3, the sign switches when the number of quarks is more than nine 
times the number of baryons. 

Another surprising result in Table [n] concerns the pK~ balance function. Even though the K~ meson has an 
anti-up quark, the QGP component is negative. This derives from \x a b being larger for us than for uu. For the range 
of parameters explored here, the hadronization component of the pK~ balance function was always positive. This 
makes it easy to recognize the existence of both the QGP and hadronization components, and if such a structure were 
observed experimentally, it would be difficult to explain without a two-component picture of quark production. 

The upper two tables in Table |III| show the dependence of the weights for variations of the baryon suppression 
Jb, which scales the final baryon yields relative to thermal yields. The values of fs roughly span the range of 
uncertainties from the experimental measurement. Whereas the default value of fs was set to 0.6, Table [TTT] shows 
results for fg = 0.5 and fg = 0.7. The number of quarks per unit rapidity in the QGP just before hadronization is also 
uncertain, hence a range of quark numbers is explored. Bracketing the default ratio of quarks before hadronization to 
final-state hadrons of 0.85, results for 7T, quar k s /«had = 0.7 and n quar k s /ft-had = 1-0 arc shown in the bottom two tables. 
The pp balance function is especially sensitive to both numbers. The hadronization peak is strengthened by raising 
fs-, or by lowering n qU arks/^had- For lower baryon yields, or higher quark densities, the hadronization peak becomes 
smaller, and can even become negative. These would lead to a dip in the pp balance function at small relative rapidity, 
which would both provide striking evidence of the two-wave nature of quark production, and suggest that the QGP 
was rather quark-rich. This latter conclusion could be better strengthened by better measurements of baryon yields, 
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which differ from coilaborations by several tens of percent. 



BLAST- WAVE PREDICTIONS FOR BALANCE FUNCTIONS 



Once one has calculated the weights described in the previous section, one can calculate the balance function between 
any two species in coordinate space given the characteristic widths of the distributions, <T(qgp) and cr^ad)- There is no 
firm understanding of the scale C(qgp), as the value depends on the microscopic details of how quark-antiquark pairs 
are created in the pre-thermalization stage. In a flux-tube picture, the quarks are pulled apart longitudinally, using 
the tube's energy to create the particles. From balance functions of pp collisions, one would estimate ct^qqp) > 0.5 
units of rapidity. Even if the balancing quark pairs are created atop one another, one would still expect the range 
to be near 0.5 once one accounted for diffusion, which spreads logarithmically with the time [JJ. The characteristic 
spread for CT(had) should be determined by the diffusion that occurs after hadronization. Although the time from 
hadronization (~ 7 fm/c) to breakup (~ 14 fm/c) is similar as the time from creation to hadronization, the diffusion 
width grows logarithmically with time, and the post hadronization diffusion should be < 0.2 units of rapidity. Since 
the spread from final-state thermal motion is likely larger than cr(h a d), the choice of 0.1 vs 0.2 for the width should 
not strongly impact the results. Since the purpose of this section is to provide an example providing a crude idea of 
what one might suspect, the values are picked with some arbitrariness to be <7(qgp) = 0.6, <7(had) = 0.2, i.e., 

(had) (QGP) 

B {Ar)) = ^1 e -(A„)V(^ ad) ) + e-^/KW. (25) 

(27r) 1 /^ C r (had) (27r) 1 /V (QGP) 

Unfortunately, the balance function is not measured in coordinate space. The mapping of r\ — > y has a spread from 
the thermal motion of the particles at breakup. For pions this can be a half unit of rapidity, whereas for protons 
the thermal spread is only a few tenths. Additionally, particles decay. To include both decays and the thermal 
spread, the correlations in coordinate space were overlaid onto a simple blast-wave parameterization. The blast-wave 
parameterization models the collective and thermal motion by assuming that the radial flow grows linearly in radius, 
ri/fmax, with u max = 1.0, and that the breakup temperature is 100 MeV. Decays of unstable particles 
(lambdas, neutral kaons, Sigmas, Cascades and Omegas) were also accounted for by a Monte-Carlo simulation. 

Three of the resulting balance functions are presented in Fig. [T] and are broken down by components. The 7r + 7r~ 
balance function is dominated by the hadronization component, with the QGP component contributing to the tail. 
The contribution from final-state decays is small, but non-negligible. Due to large thermal spread for pions, it is 
difficult to distinguish the two components. 

The pp balance function, displayed in the middle panel of Fig. [JJ is dominated by the QGP component, with the 
hadronization component being small or negative. The calculation includes the contribution to protons from weak 
decays, and the hadronization component from hyperon balance functions makes the hadronization component more 
negative. Higher quark densities in the QGP or lower final-state baryon yields push this component toward being 
negative. If it is negative, as in the case of the default calculation in Fig. [T] the resulting balance function has a 
plateau or perhaps a dip at small Ay. The existence of such a dip would provide striking evidence for the two-wave 
nature of charge production. If the hadronization component were zero, one could still see evidence of two components 
by comparing with the 7r + 7r~ balance functions. Since the pp balance function is dominated by the QGP component 
while Tr + n~ balance function is driven by the hadronization component, one could perform a single- wave fit to the 
width of the balance function in coordinate space a^. One would expect the width for pp to be significantly larger 
than that for it + tt~ . The width, for 7r + 7r - in coordinate space has been determined by a blast-wave analysis in 
[IT]. By using blast-wave parameters fit to spectra and elliptic flow observables, the analysis determined that the 
width of the balance function in coordinate space, assuming a single scale for the charge correlation, was ~ 0.22 
for the most central collisions. By performing the blast-wave fit, the contribution from final-state thermal motion 
was effectively subtracted to find the width in coordinate space. This width fell for increasing centrality from 0.6 to 
nearly 0.2. For the pp balance function one might see stay roughly constant with centrality. 

The pK~ balance function, from the lower panel of Fig. Q, offers yet more promise for demonstrating the two- 
component nature of the balance function. Since the weights of the two components have different signs with the 
stronger component having the smaller width, one finds both positive and negative regions of the balance function. If 
there were only one component, or if the two components had similar widths, this behavior would not ensue. Further, 
since protons and kaons are more massive, the thermal spread is reduced and the reduction in smearing allows more 
resolving power into the correlations in coordinate space. The negative dip for Ay ~ 1.0 might be reduced if the 
two waves of charge production are not well separated. For instance, if hadronization was more gradual the narrow 
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f B — 0-5, w quarks /n had = 0.85 





P 


A 


E+ 


E" 




s 


CT 


7T + 


K+ 


p 


0.535,-0.164 


0.586,-0.260 


0.591,-0.246 


0.581,-0.273 


0.643,-0.342 


0.638,-0.355 


0.694,-0.438 


0.005, 0.013 


-0.052, 0.096 


A 


0.221,-0.098 


0.278,-0.132 


0.278,-0.132 


0.278,-0.132 


0.334,-0.167 


0.334,-0.167 


0.391,-0.201 


-0.000, 0.000 


-0.057, 0.034 


g" 


0.089,-0.037 


0.111,-0.053 


0.113,-0.049 


0.110,-0.057 


0.135,-0.065 


0.133,-0.069 


0.156,-0.081 


0.002, 0.004 


-0.022, 0.016 


E+ 


0.088,-0.041 


0.111,-0.053 


0.110,-0.057 


0.113,-0.049 


0.133,-0.069 


0.135,-0.065 


0.156,-0.081 


-0.002,-0.004 


-0.023, 0.012 




0.056,-0.030 


0.077,-0.038 


0.077,-0.037 


0.076,-0.039 


0.098,-0.046 


0.098,-0.047 


0.119,-0.054 


0.000, 0.001 


-0.021, 0.009 




0.055,-0.031 


0.077,-0.038 


0.076,-0.039 


0.077,-0.037 


0.098,-0.047 


0.098,-0.046 


0.119,-0.054 


-0.000,-0.001 


-0.021, 0.008 


n+ 


0.011,-0.007 


0.016,-0.008 


0.016,-0.008 


0.016,-0.008 


0.021,-0.010 


0.021,-0.010 


0.027,-0.011 


-0.000,-0.000 


-0.005, 0.001 




0.121, 0.319 


-0.000, 0.000 


0.242, 0.639 


-0.242,-0.639 


0.121, 0.319 


-0.121,-0.319 


-0.000, 0.000 


0.242, 0.639 


0.121, 0.319 


K~ 


-0.248, 0.459 


-0.718, 0.437 


-0.694, 0.502 


-0.742, 0.373 


-1.164, 0.480 


-1.188, 0.415 


-1.634, 0.458 


0.024, 0.064 


0.470, 0.022 










Jb — 0.7, n quar k s /nhad — 


0.85 










P 


A 


E+ 


E" 






rr 


7T + 


K+ 


P 


0.375, 0.004 


0.411,-0.092 


0.418,-0.074 


0.405,-0.111 


0.455,-0.170 


0.448,-0.189 


0.492,-0.266 


0.007, 0.019 


-0.037, 0.096 


A 


0.155,-0.035 


0.217,-0.067 


0.217,-0.067 


0.217,-0.067 


0.279,-0.099 


0.279,-0.099 


0.341,-0.131 


-0.000, 0.000 


-0.062, 0.032 


E~ 


0.063,-0.011 


0.087,-0.027 


0.089,-0.021 


0.085,-0.032 


0.113,-0.037 


0.110,-0.042 


0.136,-0.052 


0.002, 0.006 


-0.024, 0.016 


E+ 


0.061,-0.017 


0.087,-0.027 


0.085,-0.032 


0.089,-0.021 


0.110,-0.042 


0.113,-0.037 


0.136,-0.052 


-0.002,-0.006 


-0.026, 0.010 




0.039,-0.015 


0.064,-0.023 


0.065,-0.021 


0.063,-0.024 


0.089,-0.029 


0.089,-0.031 


0.114,-0.037 


0.001, 0.002 


-0.025, 0.008 




0.039,-0.016 


0.064,-0.023 


0.063,-0.024 


0.065,-0.021 


0.089,-0.031 


0.089,-0.029 


0.114,-0.037 


-0.001,-0.002 


-0.025, 0.006 


n+ 


0.008,-0.004 


0.014,-0.005 


0.014,-0.005 


0.014,-0.005 


0.020,-0.007 


0.020,-0.007 


0.027,-0.008 


-0.000, 0.000 


-0.006, 0.001 




0.118, 0.317 


-0.000, 0.000 


0.236, 0.634 


-0.236,-0.634 


0.118, 0.317 


-0.118,-0.317 


-0.000, 0.000 


0.236, 0.634 


0.118, 0.317 


if" 


-0.126, 0.331 


-0.560, 0.291 


-0.536, 0.355 


-0.584, 0.227 


-0.971, 0.315 


-0.994, 0.251 


-1.405, 0.275 


0.024, 0.064 


0.434, 0.040 










fB = 0.6, riquarks/flhad = 


= 0.7 










P 


A 


E+ 


E" 




H~ 




7T + 


K+ 


P 


0.363, 0.012 


0.399,-0.076 


0.404,-0.059 


0.394,-0.093 


0.440,-0.147 


0.435,-0.165 


0.476,-0.236 


0.005, 0.017 


-0.036, 0.088 


A 


0.151,-0.029 


0.200,-0.052 


0.200,-0.052 


0.200,-0.052 


0.248,-0.075 


0.248,-0.075 


0.297,-0.097 


0.000,-0.000 


-0.049, 0.023 


E" 


0.061,-0.009 


0.080,-0.021 


0.081,-0.016 


0.078,-0.026 


0.100,-0.027 


0.099,-0.032 


0.119,-0.039 


0.001, 0.005 


-0.019, 0.012 


E+ 


0.060,-0.014 


0.080,-0.021 


0.078,-0.026 


0.081,-0.016 


0.099,-0.032 


0.100,-0.027 


0.119,-0.039 


-0.001,-0.005 


-0.020, 0.007 




0.038,-0.013 


0.057,-0.017 


0.057,-0.016 


0.057,-0.019 


0.076,-0.020 


0.076,-0.021 


0.095,-0.024 


0.000, 0.001 


-0.019, 0.004 




0.038,-0.014 


0.057,-0.017 


0.057,-0.019 


0.057,-0.016 


0.076,-0.021 


0.076,-0.020 


0.095,-0.024 


-0.000,-0.001 


-0.019, 0.003 


n+ 


0.007,-0.004 


0.012,-0.004 


0.012,-0.004 


0.012,-0.004 


0.017,-0.004 


0.017,-0.004 


0.022,-0.005 


-0.000,-0.000 


-0.005, 0.000 




0.098, 0.339 


0.000,-0.000 


0.196, 0.678 


-0.196,-0.678 


0.098, 0.339 


-0.098,-0.339 


-0.000,-0.000 


0.196, 0.678 


0.098, 0.339 




-0.144, 0.353 


-0.517, 0.242 


-0.497, 0.310 


-0.536, 0.173 


-0.869, 0.199 


-0.889, 0.131 


-1.241, 0.088 


0.020, 0.068 


0.372, 0.111 










fs — 0.6, riquarks/nhad - 


= 1.0 










P 


A 


E+ 


E" 




3~ 


Q- 


7T + 


K+ 


P 


0.519,-0.144 


0.570,-0.247 


0.577,-0.232 


0.563,-0.262 


0.629,-0.336 


0.622,-0.351 


0.681,-0.440 


0.007, 0.015 


-0.052, 0.104 


A 


0.215,-0.093 


0.285,-0.137 


0.285,-0.137 


0.285,-0.137 


0.355,-0.181 


0.355,-0.181 


0.425,-0.225 


0.000,-0.000 


-0.070, 0.044 


E" 


0.087,-0.035 


0.114,-0.055 


0.116,-0.050 


0.112,-0.059 


0.143,-0.070 


0.141,-0.075 


0.170,-0.090 


0.002, 0.005 


-0.027, 0.020 


E+ 


0.085,-0.040 


0.114,-0.055 


0.112,-0.059 


0.116,-0.050 


0.141,-0.075 


0.143,-0.070 


0.170,-0.090 


-0.002,-0.005 


-0.029, 0.015 




0.055,-0.029 


0.081,-0.042 


0.082,-0.040 


0.081,-0.043 


0.109,-0.053 


0.108,-0.054 


0.136,-0.065 


0.001, 0.001 


-0.027, 0.012 




0.054,-0.030 


0.081,-0.042 


0.081,-0.043 


0.082,-0.040 


0.108,-0.054 


0.109,-0.053 


0.136,-0.065 


-0.001,-0.001 


-0.028, 0.011 


n+ 


0.011,-0.007 


0.017,-0.009 


0.017,-0.009 


0.017,-0.009 


0.024,-0.012 


0.024,-0.012 


0.031,-0.014 


-0.000,-0.000 


-0.007, 0.002 




0.140, 0.297 


0.000,-0.000 


0.281, 0.594 


-0.281,-0.594 


0.140, 0.297 


-0.140,-0.297 


0.000,-0.000 


0.281, 0.594 


0.140, 0.297 




-0.206, 0.415 


-0.738, 0.463 


-0.710, 0.523 


-0.766, 0.403 


-1.241, 0.571 


-1.269, 0.512 


-1.773, 0.620 


0.028, 0.060 


0.532,-0.048 



TABLE III. Weights ™^ GP ' and w^ d ^ as in Table jll| except for a higher and lower yields of baryon number (adjusted by the 
scaling factor fs), and quark densities before hadronization, n qU arks- The upper two tables show weights for variations of /s 
from the default value of 0.6 to 0.5 and 0.7, while the lower two tables display results for varying n qU a r k s /^had from the default 
value of 0.85 to 0.7 and 1.0. The pp weights are considerably sensitive to both the baryon yield and the input baryon density. 
Either higher final-state baryon yields, or lower quark densities strengthen the hadronization peak in the pp balance function. 
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peak would have a long non-Gaussian tail which could overwhelm the smaller negative contribution from the QGP 
component. However, even in this case the overall width of the pK~ balance function would be narrowed by the second 
component. Thus, the signal for two waves of charge production which are only semi-distinct would be a narrow pK~ 
balance function, whose width might even be narrower that what one would predict from a single-wave model with 
zero width after one corrects for thermal broadening. The difficulty with pK~ balance functions comes from the 
fact that they are smaller, by nearly an order of magnitude, than the tt + tt~ balance functions, and thus require 
high-statistics data sets. Fortunately, both STAR at RHIC and ALICE at the LHC provide both high statistics, and 
due to the installation of large-coverage time-of-flight detectors, large acceptances for identified particles. 

SUMMARY 

A central feature of the canonical picture of the chemical evolution of the quark-gluon plasma is the two- wave nature 
of quark production. Investigating balance functions over a large range of species pairs provide the means to test this 
hypothesis in great detail. Once baryon production is better understood, the only parameter affecting the calculation 
of weights is the quark density in the QGP. The blast-wave parameters used to model the thermal broadening of 
the balance-function structures are already well determined by spectra. This leaves three parameters, ^ qU arks/^hadj 
cr (had) an d C(qgp) f° r fitting the entire array of balance functions. If one were to also question the assumption that 
the strange quark density was close (~ 90%) to the up or down quark density, one would add a fourth parameter. 

The existence of two waves of charge production could have several clear signatures. 

• The width of the 7r + 7r~ balance function in Ar/ (coordinate space rapidity) should be small for central collisions. 
This has been reported in [TT] . 

• In central collisions the width of the pp and K + K~ balance functions in Ar) should be larger than that of the 
7r + 7r~ balance function. Whereas the width for pions has been observed to shrink with centrality, these widths 
may well stay fixed, or even broaden for increasing centrality. 

• The pp balance function could have a plateau or even a dip at small Ay. 

• The pK~ balance function should be narrower than can be fit with a single- wave picture, and might dip negative 
for Ay ~ 1.0. 

Aside from qualitatively demonstrating the two-wave nature of quark production, the numerical parameters one 
might extract by fitting to data are also of high interest. To date, there has not been a convincing means for extracting 
the number density of quarks in the plasma, n quar k s > from experiment. Determining the width, C(qgp)j would provide 
insight into the dynamical mechanism for the creation and diffusion of quarks in the plasma. 

Potentially, the most important implication of charge balance function would be to quantitatively constrain the 
charge correlations in the QGP. For this study, the density of quarks was varied, which then determined the magnitude 
of the diagonal components of g a b(Ar]) in the QGP. Several of the hadronic balance functions were then found to be 
sensitive to this number. Additionally, there was an explicit assumption that the off-diagonal elements were zero in the 
QGP. This would not be the case if quark-antiquark pairs, such as pionic fluctuations, made significant contributions 
to the entropy of the QGP. Observing that the off-diagonal elements were small or zero, would make a strong case 
that quarks are the dominant quasi-particles in the QGP. In principle, one could extend the ideas presented here and 
vary the off-diagonal elements to determine the ranges to which they are constrained by experiment. 

The calculations presented here are somewhat schematic in nature, and can be improved during the coming years. 
Most immediately, the list of resonances considered was small, and omitted the short-lived hadronic states such as the 
p or A, which should have some measurable effect [TH]. Although almost all particles have a decay in their history, 
on the order of 10% of the charged particles produced at RHIC come from the decays of neutral resonances, other 
than the weak decays accounted for here, where both charges escape untouched from the decay. Thus, several of the 
weights might be affected at the 10% level in a more thorough calculation. If the data indeed seems addressable with 
this schematic model, one could consider more sophisticated models of quark production, diffusion, hadronization and 
emission. 
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FIG. 1. (color online) Balance functions for tv + tt~ are shown as a function of relative rapidity (black circles) in the upper 
panel as calculated with the default values given in Table |TT| The widths assumed are <7( had ) = 0.2, <7(qgp) = 0.6. The blast- 
wave model is used to map B(Arj) to B(Ay). The hadromzation component (green squares) is larger and narrower than the 
QGP component (red upward triangles), due to the fact that most quark-antiquark pairs are created at hadronization. The 
contribution from weak decays (blue downward triangles) to two pions is also shown. Since the QGP contribution is rather small 
and since the thermal spread for pions is large, one can not distinguish the two components, but must be content with the width 
being consistent with being dominated by a single narrow contribution. In contrast, the pp and pK~ balance functions shown 
in the middle and lower panels clearly illustrate the two-component nature. Because the pp balance function has contributions 
with opposite signs, one can expect to see a plateau or even a dip at small relative Ay. Since the pp balance function is 
dominated by the QGP contribution while the 7r + 7r~ balance function is largely driven by the hadronization contribution, one 
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